function ms=compute_gould_stats5m(t1,t2)
% ms=compute_gould_stats5m(t1,t2) computes stats
% for the 5-min wind stress and heat flux time
% series for the period [t1:t2].
% 4/11/01

load gould_met5m
m=find(t>=t1 & t<=t2);
d2r=pi./180;
tx=stress.*cos(d2r*wd);
ty=stress.*sin(d2r*wd);
ms=[mstats(tx(m));mstats(ty(m));mstats(qnet(m));mstats(qsw(m));mstats(qlw(m));mstats(qsen(m));mstats(qlat(m))];
